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In this paper, it is the first time to construct a complete post-Newtonian (PN) model of a rigid 
body by means of a new constraint on the mass current density and mass density. In our PN rigid 
body model most of relations, such as spin vector proportional to the angular velocity, the definition 
on the moment of inertia tensor, the key relation between the mass quadrupole moment and the 
moment of inertia tensor, rigid rotating formulae of mass quadrupole moment and the moment of 
inertia tensor, are just the extension of the main relations in Newtonian rigid body model. When 
all of 1/c 2 terms are neglected, the PN rigid body model and the corresponding formulae reduce 
to Newtonian version. The key relation is obtained in this paper for the first time, which might be 
very useful in the future application to problems in geodynamics and astronomy. 
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The idea of Newtonian rigid body has been used to treat the rotation of astronomical bodies up to now. For 
example, most Newtonian treatments of Earth's rotation are based on an accurate rigid body theory (such as SMART 
97 0,0) plus the perturbative argument from elasticity, the oceans, the atmosphere, the core and so on 0. The 
idea of a rigid body in Newtonian theory is very powerful in introducing three principal axes of a body, the spin 
axis, rotation axis and figure axis without ambiguity. Then it makes the problem much simpler since there exists the 
key simple relation between the quadrupole moment and the moment of inertia tensor in a rigid body. But even in 
Newtonian theory the concept of a rigid body is only an ideal one, because there is no real rigid body in the physical 
world ,4]. Owing to the modern high accuracy requirement, the Newtonian theory has to be replaced by Einstein's 
general relativity (at least its first post-Newtonian (PN) approximation). The problems of the post-Newtonian rigid 
body have been discussed for a long time since Born's kinematical rigidity (see Dixon's review |5|). The kinematical 
rigidity is dependent on the internal velocity distribution within the body while not considering the stress and energy 
fluxes contribution on the energy- momentum tensor T a ". Dixon 0, Thorne & Giirsel 0, Klioner 0,110 and Soffcl 
[To) have done in a much better way, so-called dynamical rigidity way, in which T af3 of the body and the gravitational 
field caused by the body satisfy some certain interdependency. The interdependency is not the same for different 
authors. However all of the different interdependencies make the PN spin S l proportional to the angular velocity £l a 
and define a relativistic moment of inertia tensor. Certainly, the interdependency for the Newtonian rigid body just 
corresponds to the simple formula between mass density and current density (shown in Eq. Q). But none of them 
obtains the key simple relation between the PN quadrupole moment M ao and the PN moment of inertia tensor I ao 
like Newtonian one as shown in Eq. @. Some scientists even assert that such a key relation is invalid in general 
relativity . Therefore the idea of PN rigidity almost has not been directly applied to practical problems up to 
now. We have a different opinion. We think because no one has discovered a suitable interdependency inside the 
energy-momentum tensor and the gravitational field before, then the key relation between PN M ao and PN I ao is not 
found. In this paper we present a suitable new interdependency to obtain the PN l ab and the key relation between PN 
M a b and PN I ab similar to Newtonian one. It is the first time to solve the rigid body problem on the post-Newtonian 
level. Recently we suggested another interdependency inside the energy-momentum tensor and gravitational field of 
the rigid body on PN level by means of a special gauge condition , but the special gauge condition is more or less 
speculative and difficult to be commonly accepted. In this paper we totally discard the special gauge condition, and 
present a general expression of the rigid spin with or without external field (free precession). 

First let us recall the basic aspects of the Newtonian rigid body. We take E and E a (= T,V a ) as the mass density and 
the mass current density of a rigid body A respectively. Then the mass multipole moments Mj, and spin S a of body 
A are defined as M L = f A d 3 XZX L and S a = e abc J A d 3 XX b Z c , where X L is an abbreviation of X Kil X 12 ■ ■ ■ X ll> , 
in which ij (j — 1, 2, • • ■ , I) is a spatial index, the angular brackets mean "symmetrize and take the trace-free part" 
(STF). In Newtonian rigid body the rotational angular velocity fi is independent of spatial coordinate, we have 

E a = e abc Zn b X c . (I) 

Insert Eq. (JTJ into the expression of spin, we have the linear relation between the spin and angular velocity S a — I a b^ b , 
where the moment of inertia tensor I a f, is I a b = ha = Ja d 3 XYj(X 2 8 ao — X a X b ). 
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The mass quadrupole moment and the moment of inertia tensor satisfy the key relation 

M a b = —I a b + ^ablcc ■ 



(2) 



By means of the continuity equation <9t£ + d a Y, a — 0, the time derivative of the moment of inertia tensor I a b is 
proportional to the angular velocity $7 



lab — (] j r ] n ^ a ^ (^apqlqb ~l~ ^bpq-^aq)^^ • 



(3) 



M a b satisfies a similar relation as ©. Therefore, I ao and M ao like constant tensors, will rigidly rotate in space with 
the angular velocity f2. 

When we discuss PN rigidity, we will use the symbols, notations and conventions following the 1PN theoretical 
framework presented by Damour, Soffel and Xu (cited below as the DSX scheme 0,EJE1) since the DSX scheme 
is not only rather simple and complete but also describes the 1PN definition of spin in a satisfactory manner. In 
the DSX scheme a complete 1PN general relativistic celestial mechanics for N arbitrarily composed and shaped, 
rotating deformable bodies is described. Here we will briefly summarize the notation and definition in the DSX 
scheme. In the post Newtonian expansion we will always abbreviate the order symbol 0(c~ n ) simply as 0(n). A 
spatial multi-index containing I indices is simply denoted by L (and K for k indices, etc), i.e. L = A 
multisummation is always understood for repeated multi-indices SjJTl = X)j ' ' ^his— ii^im— k- Given a 
spatial vector, n 1 , its Ith tensorial power is denoted by n L = ri Ll n 12 ■■ - n H . Also, 8l = d^d^ ■ ■ ■ d^. Besides angular 
brackets the symmetric and trace-free (STF) part of a spatial tensor will be denoted by a caret when no ambiguity 



arises: STF, 



of the Cartesian metric Si 



■■ti> — 
diag( 



T^. The spatial indices, i,j = 1,2,3 are freely raised or lowered by means 
fl, +1, +1) in Cartesian coordinates. The metric is presented by means of 



potential W and vector potential W a (see Eq. (4.1) of Ref. 0). W and W a can be separated into self-part (with a 
"+") and external-part (with a bar), i.e. W = W + + W and W a = W+ +W a . The self-part W+ and W+ will be 
solved by the gravitational mass density S and the mass current density E a as: S = (T 00 + T ss )/c 2 , and S a = T 0a /c, 
through the 1PN Einstein field equation and the coordinate conditions (gauge conditions) (see Eq. (4.3) of Ref. [12j), 
where T a/3 is the energy-momentum tensor. W + and will be expanded by STF BD mass moments Ml and STF 
spin moments Sl (see Eq. (6.11) of Ref. 12]). The external part W and W can be expanded by gravito-electric 
tidal moments Gl and gravito-magnetic tidal moments Hl (se e Eq. (4.15) of Ref. [Il|). Gl and Hl are also STF 
spatial tensors dependent on time only. BD mass moments |l5j are widely accepted as the best PN mass moments, 
which have the form 



m£(t) 



4(2/ + 1) d 



d 3 XX L X 2 Y; 



d 3 XX aL Y; a 



(I > 0) . 



(4) 



(I + 1)(2Z + 3)c 2 dT 

The PN spin moment has been discussed for a long time [l4lHr|. In the DSX scheme, the expression of 1PN spin of 
body A is 



A, pn _ 



Cafcc / d XX 
A 



E c (l 



^ W A ) _^ {W +A + l dcdTZ + ) 



■-T- 



1 



2/ + 3 



Ht L N£ 



I 



l + l 



• 2c 2^Z!(/ + 2)(2Z + 5) 



(/ + 10)N b A L G? L 



+ 8(21 + 3)P b A L G± -(1 + 2)N A bL> G A L + 0(4) , 



(5) 



where Z\ = G f A d 3 X'E(T A , X')|X - X'|, W+ A = G J A d 3 X'Z c (T A , X')/|X - X'|, dot means time derivative d T , 
N L and P L are defined as (see Eq. (2.10) of [3) N L = J A d 3 XX 2 X L E and P L = J A d 3 XX aL T, a respectively. Later 
we omit to indicate the body label, A, on all quantities. In Ref. ^ij; Nl and Pl are called "bad moments". In such 
a definition of 1PN spin vector, S 1 ™ satisfies the rotational equation of motion (see Eq. (3.11) of Ref. 01)- We also 
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have 1PN continuity equation for (E, E a ) (see Eq. (5.6b) of Ref. [12() 

1 



d T E + d a E a = ^(d T T bb - Ed T W) + 0(4) . 



(6) 



Those are the equations which are taken from the DSX scheme and will be used in the following discussion on the 
PN rigidity. 

The definition of the 1PN rigid body has to agree with the Newtonian rigid body when 1/c 2 terms are neglected. 
In the 1PN rigid body the angular velocity should be independent of the local coordinate X a of body A. In the 
DSX scheme we substitute E and E c for the energy-momentum tensor T af3 , therefore the interdependency inside the 
energy-momentum tensor and gravitational field in Ref. 00)13 might be replaced by inside E, E c and gravitational 
field. We expect that the interdependency will produce equations similar to and J3J) on the PN level. On the PN 
level it is sufficient to replace T Q/3 by E, E c and their derivative 01- Before a further discussion on the rigid body, 
the rigid BD moments and the rigid 1PN spin vector should be considered. In the rigid body, since Pl and Nl are 
PN terms we can substitute the Newtonian relations (Eq.JQ) and Newtonian continuity equation) for the definitions 
of Nl and Pl . It is easy to prove 



P. 
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<L> 



21 + 1 



-N, 



<L> 



Lemma 1 The rigid BD mass moments of rigid body A (Eq. Q) can be simplified as 



Mr 



d 3 XX <L> 



1 + 9 



c 2 V2(Z + 1)(2Z + 3) 



X 2 E 



•0(4). 



(7) 



(8) 



Proof of Lemma 1. Starting from Eq. (the definition of BD mass moment) and replacing the third term in 
the RHS of Eq. J3J by Eq. J7J). The prove will be carried out directly. In fact only the relativistic quadrupole 
moment is interested in, because in the solar system all of the relativistic higher multipole moments are too small 
to be considered in any modern measurement in the anticipation of the day. Therefore we only take the relativistic 
quadrupole moments into account: 



M ab = I d 3 XX <ah> 

I A 
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42c 2 



-X^ 



We define the PN part of the mass density E 



pn = ||X 2 S and the total mass density 
— _ , Epn 



(9) 



(10) 



Lemma 2 The rigid PN spin vector of body A (see Eq. 0) can be reduced to 



;>PN 



e abc I <PXX b { E 

E 
7! 



E 



-e cde n d d e z+ + -e edf n d xfd ce z^ 



r 2 



Ae cde n d X e X <L> G L (T) 



1 



Z + 10 



2(Z + 2)(2Z + 5) 



X L X 2 G cL + 



1 + 2 
I + 10 



iX 



<dL> 



2(l + 2)(2l + 5) 



d T (lnZ)X L X 2 G cL 



(11) 



Substituting Eqs. and Q for the PN part of Eq. JSJ, integrating by parts, using the surface integration for whole 
body A to be zero and taking some STF formulae, we obtained Eq. ((TIT) . 
We define the PN self-part and PN external part of current density as 



E s c elf = E (le cde Q, d d e Z + + ^e edf n d Xfd ce zA , 



y- 

l>0 



4e cde n d X e X <L> G L (T) 
I + 10 



1 



2(Z + 2)(2Z + 5) 



X^X z G cL 



1 + 2 
I + 10 



£cedX <dL> H e L 



2(Z + 2)(2Z + 5) 



<9 T (lnE)X L X 2 G cL 



(12) 



(13) 
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Both of Sggif and E£ xt as well as E c itself are spatially compact-supported. When tidal moments Gl and Hl are 
equal to zero, E c and Eg elf /c 2 form the PN self-part of spin vector. We can define 

S c = E c + + ^ft , (14) 



then Eq. (|1 If) becomes 

S a PN = e abc I d 3 XX b ^ . ( 1.5 ) 



By comparing the Newtonian definition of spin with Eq. (fT3|) . £ c is the fully PN quantity. 
We add -4 [<9t(E pn ) + <9 Q (E£ clf + EJJ] to both sides of Eq. © and have 



^ „ 1 



a T s + a a s = - [d T T bb - m T w + 9 T Sp N + d a (^ cU + e^j] . (ie) 



Now we construct the model of the PN rigid body by constraint on E and S to satisfy 
1 



2c 2 



X a [d T T bb - Hd T W + d T E PN + d a (E s a olf + E£J] = e afcc ft b X c E + 0(4) . (17) 



The relation (Eq.JTJJ) is our most important assumption for PN rigid body. Considering that £ c and E are expressed 
by E c and E, which are related to T a @ in the DSX scheme, then Eq. 117(1 is also constraint to T a @ . When 1/c 2 terms 
are neglected, Eq. 1|17[) goes to Eq. QJ. Later we will see that only in this model PN mass quadrupole moments 
and the moment of inertia tensors satisfy the similar Newtonian key relations as in Eq. J2J). We were not surprised 
with the appearance of the time derivative of E in Eq. I|17l) . since in the DSX scheme T a @ can be fully represented 
by E and E a and their derivatives without difficulty [17j. In the interdependencies described by Thorne and Giirsel 
and Klioner ||, they have their own models of the rigid body by another constraint to T a @. By comparing the 
constrained equations (to see Eq.(A7) in 6] (or Eq.(7) in ||) and Eq.(8) in @) with Eq. JlZj), we see that Eq. fT7| 
is more complicated, but still reasonable. 

Substituting Eq. (HI for Eq. (JUJ), we obtain the linear relation between the PN spin vector of the rigid body and 
the angular velocity: 



S™ = I ab n b + 0(4) , (18) 



where the moment of inertia tensor is 



I ab = ha= I d 3 X(S ab X 2 - X a X b )H + 0(4) , (19) 

J A 



in which E is defined in Eq. HlOfl . 

By comparing Eq. I|19fl with Eq. 10, we have 

M ab = -I ab + ^5 ab I cc + 0(4) . (20) 

Eq. (|20|l is the key relation between the PN mass quadrupole moment (rigid BD moment) and the PN moment of 
inertia tensor. It is just this relation that makes the model of the rigid body very useful and applicable on the PN 
level as shown in the Newtonian case. This is the first time that we have obtained the PN key relation in this paper. 
Making use of the extended 1PN continuity equation Eq. (|16f) . we immediately have 

i ab = -?=I ab = (e apq I qb + e bpq I aq )n p + 0(4) . (21) 

PN M ab satisfies a similar relation as l(21|l . From Eq. H21[) the behavior of the PN I ab (and also PN M ab ) in our rigid 
model are just like the Newtonian version (Eq. ©), i.e. I ab and M ab rigidly rotate as a whole. 

Eq. (|21|l means that we always can introduce a rotation matrix Pi a (T), which is a time-dependent orthogonal matrix 
and transform the PN reference system (RS) to the corotating reference system with rigid body (RS + ). Pi a (T) can be 
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constructed by the rotational angular velocity fl a of rigid body according to the relation: Sl a (T) = | 
(H3)- I n the new corotating coordinates we get 



€abcPib(T)Pic(T) 



(22) 



where ly = P ia P bIab- Pia satisfies the following relations: P ia Pj a = PaPb = <5 a b and dP ta /dT= e abc VL b P lc (here 
we use fr to substitute for a;- 7 in [13). The proof is easy by means of i|21[l. Eq. 122(1 shows it is possible to introduce 
the PN Tisserand reference system. 

At last, we should emphasize that the calculation of the PN moment of inertia tensor Eq. (|19() is not too difficult, 
although the constraint relation on E and S in Eq. (|17|l in the model of the PN rigid body is complicated. In practical 
problems, from our PN rigid spin, the PN moment of inertia tensor (Eqs. i|18|l - (123)) it is possible to define the three 
principal axes of the body, the spin axis, rotation axis and figure axis as described in the Newtonian theory, which 
we will discuss in a separate paper in the future. 

In conclusion, the rigid BD (PN) mass multipole moments Eq. {HJ and the rigid PN spin moment Eq. I|15(l are 
discussed in this paper. We successfully have constructed a new PN model of the rigid body in which the constraint 
on S and £ satisfies Eq. (II Til . Our PN rigid body model will reduce to the Newtonian one when all of 1/c 2 terms are 
neglected. Most of relations in our PN rigid body model, such as the spin vector proportional to the angular velocity 
f2 (Eq. (|18|l ). the definition on the moment of inertia tensor (Eq. (|19(0 . the key relation between the mass quadrupole 
moment and the moment of inertia tensor (Eq. J20I)), the rigidly rotating formulae of I a b and M a b (to see Eq. (I21f) ') 
are similar to the Newtonian rigid body model where the corresponding relations are mentioned at the beginning 
of this paper. Especially, the PN key relation between M Q f, and I a b might be applied to the practical problems in 
geodynamics and astronomy in the future, e.g. the discussion on the relativistic effects of the nutation and precession. 
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